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We consider proximity effect coupling in Superconducting/ Antiferromagnetic/Superconducting (S- 
A-S) sandwiches using the recently developed SO (5) effective theory of high temperature supercon- 
ductivity. We find that, for narrow junctions, the A region acts like a strong superconductor, and 
that there is a critical junction thickness which depends on the effective SO(5) coupling constants 
and on the phase difference across the junction, at which the A region undergoes a Freedericksz-like 
transition to a state which is intermediate between superconductor and antiferromagnet. For thick 
junctions, the current-phase relation is sinusoidal, as in standard S-N-S and S-I-S junctions, but for 
thin junctions it shows a sharp break in slope at the Freedericksz point. 



Zhang has recently developed a theory 0j which unifies 
d-wave superconductivity (S) and antiferromagnetism 
(A) on the basis of an underlying SO (5) symmetry. 
The S and A order parameters are combined into a 5- 
dimensional superspin, and the high energy physics of 
these superspins is postulated to be rotationally symmet- 
ric. At low energies this 5*0(5) symmetry is broken by 
a chemical-potential-dependent anisotropy which favors 
the A state for /i < /i c or the S state for \l > \l c . This 
implies that, at low temperature, there is a "soft direc- 
tion" for perturbations of a stable d-wave superconductor 
toward antiferromagnetism. Similarly, the appropriate 
perturbation, applied to a stable A material, will tend to 
drive it into the S state. By analogy to the proximity ef- 
fect in conventional superconductors, it is clear that the 
relevant perturbing field is provided by proximity of an A 
material to an S material. Moreover, in a sandwich S-A-S 
configuration, this proximity effect would be expected to 
provide a mechanism for Josephson coupling of the two 
S regions. We also note that one approach to practical 
high-Tc Josephson junctions involves the use of barriers 
made from the cuprates near the S/ A transition. 
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FIG. 1. Geometry of the suggested junction 

In this paper we present analytic and numerical re- 
sults for the properties of the S-A-S Josephson junction 
system, shown on Fig.|], in terms of SO (5) continuum 
theory in which the spatial variation of the order pa- 
rameter is one dimensional. We obtain analytical results 
for the critical current as a function of thickness and nu- 



merical results for the current-phase relation for different 
thicknesses. 

We find that, when the S layers are strongly super- 
conducting, thin A layers are driven completely super- 
conducting by the field of the adjacent S layers, and the 
SO (5) order parameter lies completely in the supercon- 
ducting plane. Beyond a critical barrier thickness, we 
find that the order parameter in the junction starts to 
tip back toward the antiferromagnetic plane, in a fashion 
precisely analogous to the Freedericksz transition in liq- 
uid crystals. Twisting the superconducting phase, which 
causes a current to flow through the junction, is anal- 
ogous to twisting the nematic director at the walls. A 
sufficiently large twist will drive the system through the 
Freedericksz transition resulting in a distinctive, non- 
sinusoidal current-phase relation for an S-A-S junction. 

Our results clearly demonstrate that, within 50(5) 
theory, the details of Josephson coupling through an A 
barrier are qualitatively different from those of proximity 
effect junctions with conventional barriers. Hence study 
of S-A-S junctions provides a critical test of 50(5) theory. 
By the same token our calculations provide a new ba- 
sis for the interpretation of real high-Tc Josephson junc- 
tions, currently being fabricated and studied [^]-|}. 

In the spirit of 50(5) we describe the system by 
a three-component order parameter n = { 



:}, 

where the first two components are the real and imag- 
inary parts of the superconducting order parameter and 
the third component represents the antiferromagnetic 
Neel vector ( See Figure || ) . For simplicity we treat the 
Neel vector as a single component. However this com- 
ponent may be viewed as the amplitude of a spatially 
uniform 3D vector. 

According to jlj the system is described by a functional 



£(n) = £ {d,n a f 
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with the constraint n 2 — 1. As in || we assume that 
the gradient term is 50(5) symmetric. The anisotropy 
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term g is positive in the A region (so that it would be 
antiferromagnetic in the absence of proximity effects) and 
negative in the superconductor. 
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FIG. 2. SO(5) order parameter 

The superspin constraint is most naturally imple- 
mented in polar coordinates n x = cos 9 cos 4>, n y = 
cos#sin</>, and n z = sm.9. 

£(0, <f>) = f i(d^) 2 + cos 2 9 (d^) 2 } - g sin 2 9 (2) 

In all of our calculations we will assume rigid supercon- 
ducting boundary conditions n z \o — 0, and n z \d = 0. 
Strictly speaking this is only true in the case of "strong" 
superconductors and "weak" antiferromagnets: \gs\ 3> 
\gA\- However analysis of the general case shows that 
relaxing this condition does not change the qualitative 
picture. 

At this point one can easily specify the analogy be- 
tween our problem and the problem of a liquid crystal 
in a slab with anchoring walls, in an electric field. If the 
electric field is perpendicular to the walls, it will try to 
align the director of the liquid crystal along the field. At 
small voltages the field is unable to overcome the effect 
of surface pinning, and the equilibrium configuration re- 
mains uniform. However with increasing voltage the sys- 
tem will undergo a Freedericksz transition, in which the 
director begins to align along the field. More interest- 
ingly this transition is known to depend on the applied 
boundary conditions, i.e. on the relative twist of the an- 
choring directions on the two sides of the slab (the twisted 
nematic transition) |Q. 

We now show that similar effects arise in S-A-S sand- 
wiches within 50(5) theory. The role of the voltage is 
played by d^f g A j P, and the superconducting phase differ- 
ence across the junction corresponds to the twist angle 
imposed by the two anchoring walls. The S-A-S sand- 
wiches will undergo a phase transition in which the A 
region, between the two superconductors, goes from be- 
ing purely superconducting (by virtue of the proximity 
effect) into a mixed S/A state. We also show that, suffi- 
ciently close to such a Freedericksz transition, the system 
possesses non-trivial current-phase characteristics, as a 
consequence of the transition. 

In the A region the Euler-Lagrange equations for the 
functional (0) are 
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p~r^ + p cos # sin 6* I I + 2gAsin9 cos9 — (3) 
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The boundary conditions for these equations are given 
by 



6{x = 0) = 4>(x = 0) = 
9{x = d)=0 (j)(x = d) = A$ 



(5) 
(6) 



where A$ is the phase difference between two supercon- 
ductors. 

Equation (|J) is nothing but the conservation of current. 
I s = n\d x ri2 — n%d x rii = cos 2 9-^- 



dx 



So we can write (H) as 



d 2 9 I 2 
p-nr+psm9 — s — + 2g A sm9cos9 = (7) 
dx 2 cos 4 9 

The last equation can be easily integrated once giving 



dx J cos 2 9 cos 2 9q 



with the characteristic length 

U = V p/2g A 



+ sin 2 6» (8) 
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In writing (g) we expressed the constant of integration 
in terms of the maximal value 9q that will be reached at 
x = d/2 (where ddjdx = 0). This immediately results in 
an equation for 9q 
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where uj s = I s Ca, the parameter k is defined by 



, 2 sin 2 9 cos 2 9 
uj 2 + cos 2 9q 



(10) 
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and K is the complete elliptic integral of the first kind. 
Equation ( ]To| ) should be supplemented by an equation 
for the current uj s in terms of the phase difference across 
the junction 



A$ = 2I S 
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here Hi(n, k) is a complete elliptic integral of the third 
kind. 

One can easily see that Eq. ( |l0| ) has a solution only 
when d/^A > ir/yl + wf. For smaller d the only so- 
lution will be 9q = 0, which means that the A re- 
gion remains uniformly superconducting. Even though 
antifcrromagnetism would be favored in a bulk mate- 
rial of this kind, proximity to a "strong" superconduc- 
tor forces it to be uniformly superconducting. When 
d c = tt£a/ \A + w si a second order transition occurs at 
which 6*o starts to increase as \fd— d c , so that the A 
region exhibits both kinds of order: superconductivity 
and antifcrromagnetism. It is interesting to note that a 
non-zero u> s decreases the critical width of the A region. 
This can be understood as the result of having an extra 
"torque" in the x-y plane. This result raises the very in- 
teresting possibility of choosing a width of the A region 
below the critical value at zero current d C Q = 7r£ and 
then tuning the system through the transition by simply 
passing a current through the junction! 



periodic motion in 50(5) space, permitting exploration 
of the low g-vector dynamics of 50(5) theory. 




FIG. 3. 6»o vs A<f> and L vs A(j> for an S-A-S junction. The 
onset of 8o and the discontinuity in slope of I 3 both occur at 
the Freedericksz transition. 

In Figure ^ we present such an example, for the case 
d = 0.85d c o. This figure shows that the system undergoes 
a transition when A<f> = 1.7. Below the transition 6q is 
identically zero and / s is a linear function of A<I>, as one 
would expect for a uniform superconductor. However 
above the Freedericksz transition, 6$ starts to grow and 
I s vs A<f> develops curvature. Eventually, at A<f> = ir, 
6q = tt/2 and I s goes to zero. We note that further 
interesting differences with the conventional proximity 
effect can be expected in the dynamical state at finite 
voltages. In the presence of a finite voltage across the 
junction, the full 50(5) order parameter will undergo 




FIG. 4. Go vs A<E> for junctions with different d/d c o- Notice 
that for A"I> = 7r we always have 9q = n/2. 

Figure [| shows that this feature, 8q = tt/2 when 
A$ = 7T, occurs for all widths of the A region. It may be 
understood as follows: The energy required to twist the 
superconducting order parameter by tt without changing 
its magnitude is the same as the energy required to ro- 
tate the superspin into the antiferromagnetic plane and 
back into the superconducting plane. However rotating 
the superspin into the antiferromagnetic direction allows 
the system to lower its energy because of the g-term. 




FIG. 5. Current-phase characteristics of junctions with 
different d/d C Q. 
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This effect is an interesting 5*0(5) analogue of the re- 
sult of Krotov et. al. M that superconductivity between 
antiferromagnetic stripes is suppressed for nontopological 
stripes and enhanced for topological stripes. 

Figure |B| illustrates the non-trivial current-phase char- 
acteristics of S-A-S junctions with increasing width of 
the A layer. When d < d c o they show a transition from 
linear dependence below the transition to sin-like depen- 
dence above it. Some asymmetry persists in the curves 
for d > d co , and for d ^S> d C Q they show the usual sin(A<i>) 
dependences of SIS junctions. 

It is easy to calculate the critical current of our junc- 
tions. For a given d, Eq. ([n]) does not have any solution 
for currents that are too large. The first solution appears 
at a point that corresponds to the maximum of k 2 in Eq. 
(0) §. This 

"'■max 

is given by 

Using the asymptotic forms of the elliptic functions, we 
find for Eq. @ d/(2£ A ) = ln(4/y^) which gives the 
critical current 

I. = e~ d/ ^. (13) 

So £a represents a new correlation length for supercon- 
ducting proximity effects across antiferromagnets: Ac- 
cording to H gA — 2x(m? — /^ 2 )i where \i is the chemical 
potential, and \x c is the critical value of the chemical po- 
tential at which the first order transition between the 
superconducting and antiferromagnetic states occurs. In 
deriving Eq. (^), we have assumed an 50(5) symmet- 
ric susceptibility x f° r the charge, spin and tt operators. 
According to Eq. when [i is close to fj, c (and hence 
gA is small), £4 will be large. This provides a new and 
natural explanation of the long range proximity effect 
sometimes observed in PBCO (^^]. We note, however, 
that asymmetry in the x's will generate a cut-off for £4 
in Eq. (|). For \c > Xtt we find £ max = y/p/rj where 

V = 2 f j2 c(Xc ~ Xtt). 
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